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ABSTRACT 

We present complimentary techniques to find emission-line targets and measure their properties in 
a semi-automated fashion from grism observations obtained with the Advanced Camera for Surveys 
aboard the Hubble Space Telescope {HST). The first technique is to find all likely sources in a direct 
image, extract their spectra and search them for emission lines. The second method is to look for 
emission-line sources as compact structures in an unsharp masked version of the grism image. Using 
these methods we identify 46 emission-line targets in the Hubble Deep Field North using a modest 
(three orbit) expenditure of HST observing time. Grism spectroscopy is a powerful tool for efficiently 
identifying interesting low luminosity, moderate redshift emission-line field galaxies. The sources found 
here have a median i band (F775W) flux 1.5 mag fainter than the spectroscopic redshift catalog of 
Cohen et al. (2000). They have redshift z < 1.42, high equivalent widths (typically EW > lOOA), 
and are usually less luminous than the characteristic luminosity {L'^) at the same redshift. The chief 
obstacle in interpreting the results is line identification, since the majority of sources have a single 
emission line and the spectral resolution is low. Photometric redshifts are useful for providing a first 
guess redshift. However, even at the depth of the state-of-the-art ground-based and HST data used 
here, photometric errors can result in uncertainties in line identifications, especially for sources with i 
magnitudes fainter than 24.5 ABmag. Reliable line identification for the faintest emission-line galaxies 
requires additional ground-based spectroscopy for confirmation. Of particular concern are the faint 
high EW [O II] emitters which could represent a strongly evolving galaxy population if the possibility 
that they are mis-identified lower redshift interlopers can be ruled out. 

Subject headings: Surveys; techniques: spectroscopic; methods: data analysis; galaxies: distances and 
redshifts; galaxies: high-redshift 
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1. INTRODUCTION 

One of the perennial problems in modern astrophysics 
is measuring spectral information, especially redshifts, of 
the most distant and hence, faintest sources in the uni- 
verse. With the Hubble Space Telescope (HST) and the 
Advanced Camera for Surveys' (ACS) Wide Field Cam- 
era (WFC) astronomers can now obtain reliable broad 
band photometry down to ABmag ~ 27 across the op- 
tical portion of the spectr um with a modest expendi- 
ture of telescope time (e.g. iBemtez et al]l2004f ). The in- 
strument is able to image even deeper as demonstrated 
by th e Hubble Ultra Deep Field Beckwi th et al. (HUDF 
120061 ) which has a S/N = 10 limit of ^29.9 ABmag 
for point sources with the F775W filter (|B ouwens et al. 
[2b06); in the same field and using the same filter the 
S/N = 10 detection limit for extended sources within an 
aperture having diameter of 0.4'' and 0.5'' is 29.4 and 29.0 
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ABmag, respectively (^Beckwit h et al.l 120061 : ICoe et ahl 
l2006f ). However, our limit for obtaining reliable spec- 
troscopy is much brighter. Spectra to ABmag ~ 24 are 
difficult to obt ain even with the largest ground-based 
telescopes fe.^. [Abraham et alll2004l : ICowie. et al.l[2004l : 
ICohen et al.ll2000r ). At issue is the domination of the sky 
over the signal from astronomical sources observed from 
the ground at these faint magnitudes. The background 
seen by HST is orders of magnitude weaker. However, 
slit survey spectroscopy with the now defunct STIS and 
FOS spectrographs was not feasible due to the telescope's 
small aperture and minuscule projected slit widths. Al- 
ternatively, multi-object slit less spectroscopy is and has 
been available with HST using a variety of instruments 
(FOC, STIS, NICMOS). Ah three cameras of ACS also 
have dispersing elements. The combination of the ACS 
G800L grism and the WFC is particularly noteworthy 
since it provides the widest field and highest throughput 
of all the slit less options on HST, allowing deep low- 
resolution spectroscopy with modest expenditure of tele- 
scope time. For example the three orbit integration of the 
Hubble Deep Field North (HDFN) that we discuss here 
provide S/N ~ 10 spectra in the continuum of sources 
having ABmag = 25.1 at a mean wavelength A ~ 8000 A. 

While this major new capability is welcomed, grism 
data are difficult to work with. The broad spectral cov- 
erage of G800L (X = 5600 - 9900A at 25% of peak spec- 
tral throughput, Walsh & Pirzkal 2004) results in a back- 
ground count rate that is high compared to other WFC 
filters, and indicates that wavelength variations in the 
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flat-field are a major concern. While the grism transmits 
most of its light in the first order, light from other orders 
of bright sources can also be seen. The orientation of the 
G8QQ L filter and the stro ng geometric distortion of the 
ACS (Meurer et al. 2002) results in spectral traces that 
are skewed relative to the CCD pixel grid at an angle that 
varies across the field and a wavelength calibration that 
also is field dependent. The spatially varying distortion 
makes combining dithered G800L images difficult, espe- 
cially if there are large offsets or roll-angle variations. 
Fortunately astronomers at the Space Telescope Euro- 
pean Coordinating Facility (S T-ECF) have provided cal- 
ibrat i ons of the G800L grism (|Pasquali. Pirzkal fc WalshI 
120031 : IWalsh fc Pirzkalll2004r ) as well as the aXe software 
package for e xtractin g and pro cessing slitless spectra 
(iPir zkaL Pasquali. fc Demleitnerl [2001) which relieves 
many of the pains of dealing with grism data. 

The aim of this paper is to assess techniques to process 
grism data and identify Emission Line Galaxies (ELGs). 
ELGs are particularly interesting since they mark the 
location of AGN or intense star formation and hence, 
cosmic evolution. For example, some of the most dis- 
tant galaxie s in the universe are Lya emission sources 
at 2; ^ 6.6 (iKodaira et~all l2003l : iTaniguchi erall l2005l : 
iKashikawa et al.""2006). In addition, sharp narrow emis- 
sion features should be one of the easiest spectral sig- 
natures to find in grism data. Since this is primarily a 
techniques paper, our science analysis is relatively light, 
and includes basic comparisons of the key measurable 
properties of the ELGs (luminosity, equivalent width, 
color) with other galaxy samples at a similar redshift. 
This analysis is sufficient to show that the grism pro- 
vides an efficient means for selecting statistically sig- 
nifica nt samples out to z ~ 1.5. This study is simi- 
lar to iPirzkal et al.l (|2004l ) which provides a detailed de- 
scription of how the GRism ACS Program for Extra- 
galactic Science (GRAPES) collaboration have reduced 
their G800L o bservations of th e HUDF and extracted 
source spectra. IXu et all (|2006f ) present their technique 
for identifying ELGs and th e resultant catalog of sources 
found, while Pirzkal e t al.l (j2006) discuss the morphol- 
ogy of these ELGs. Here we discuss WFC grism process- 
ing methods and tools that were developed by the ACS 
Science Team largely independently from the ST-ECF 
and GRAPES efforts and optimized for the discovery of 
ELGs. We apply these methods to observations of the 
HDFN. A brief summary of this work was presented at 
the 2005 HST Calibration Workshop (Meurer 2006 ). 

In Sec. [2] we compare our data processing to that of the 
GRAPES collaboration and describe in detail our data 
and its processing. In Sec. [3] we present two methods 
for finding ELGs, as well as our methods for assigning 
line identifications. Sec. [H presents our results includ- 
ing a list of all ELGs found, and an initial assessment of 
the statistical properties of the galaxies found. In Sec. [5] 
we compare our redshifts with other observations of the 
HDFN and determine the redshift accuracy of the grism. 
Finally, in Sec. [6] we summarize our results, discuss the 
benefits of using the grism and the additional require- 
ments for obtaining useful redshifts of ELGs. 

2. METHODS AND DATA 

Here we provide an overview of our image processing 
and object extraction and note how it differs from that of 



the GRAPES cohaboration (|Pirzkal et al.ll2004r ). Follow- 
ing that, we describe the processing of the HDFN data 
in detail. 

2.1. Comparison of image processing techniques 

The GRAPES team do minimum processing of their 
images before extracting spectra. Like us, they rely o n 
the STScI calibration pipeline CALACS (jHackl Il999f ). 
as implemented by the STScI archive to do most of 
the basic CCD processing consisting of overscan sub- 
traction, bias subtraction, dark subtraction, and gain 
correction. These steps are performed using the best 
available reference files as implemented by the STScI 
Archive. The flat-field employed by CALACS for G800L 
images is a unity flat, so in effect no pixel-to-pixel 
flat-fielding is done. The GRAPES team subtract a 
scaled super-sky frame from each image to remove the 
sky background (where the scaling is to object free re- 
gions of the image). The G800L image shifts are deter- 
mined from MULTIDRIZZLE-piocessed short exposures 
obtained through a broad-band filter at the start of each 
orbit. The MULTIDRIZZLE task is also used to produce 
geometrically-corrected G800L images, but they are used 
only to identify the cosmic rays, not for spectral extrac- 
tions. Instead, the extractions are performed by aXe on 
the individual frames after sky subtraction and masking 
of the cosmic rays. The spectra are co- added and the 
flat-fielding is performed at this stage by calculating the 
effective wavelength of the light falling on each pixel and 
interpolating between a series of broad- and narrow-band 
flats (Walsh & Pirzkal 2004) to determine the flat field 
correction most appropriate for that pixel. 

Our approach differs in a few key ways. First, we flat- 
field our G800L images with the flight flat-field for the 
F814W filter. We then process the data with the GTO 
science pipeline Apsis (Blakeslee et al. 2003) to make 
cosmic-ray rejected, aligned, combined, and geometri- 
cally corrected images. The result is one final G800L 
image which we use to extract spectra and all spectro- 
scopic quantities. There are several advantages of this 
approach. Application of the F814W flat cosmetically 
improves the images by largely removing most small- 
scale CCD blemishes. In addition it produces flatter im- 
ages, reducing the rms amplitude of the sky background 
over large scales (> 75 pixels) by a factor of two, as de- 
termined from application of the flat-field to the super- 
sky frames used by iPirzkal et al.l (|2004l ) . Since we forgo 
the use of aXe's A dependent flat-field fit, the flux scale 
varies throughout the fi eld of our images by up to ~10% 
(jWalsh fc Pirzk al"2004i). Flat-field images from WFC 
show numerous dark blemishes which are more apparent 
with decreasing A (iBohlin et al.ll2001l ). Because of the 
broad spectral response of the grism, blemishes are inac- 
curately removed with the F814W flat. For example, if 
the blue end of the spectrum of a compact source falls on 
a blemish, it will not be completely removed by flat field- 
ing and result in a spurious absorption feature. However, 
this is not a major concern here since we are concerned 
with emission lines rather than absorption features. This 
is also less of a problem for extended sources (^10 pix- 
els): since many wavelengths contribute to each pixel 
in the spectrum one can no longer assume that a single 
wavelength dominates, and hence our flat-fielding tech- 
nique should be sufficient in these instances. Using small 
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dithers can also mitigate against this happening. Com- 
bining geometrically corrected dithered data can work as 
long as there are no roll angle variations and the dithers 
are all within ~ . Then, the A scales of the first or- 
der spectra will align to within 0.025'' (~ 0.5 WFC pix- 
els) across the WFC field. Geometric correction has the 
advantage that it removes much of (but not all of) the 
spatial variation in the A calibration with the dispersion 
remaining nearly constant within each spectrum. A ma- 
jor advantage of our approach is that the geometrically 
corrected spectra are nearly horizontal. Over a spectral 
length of 75 pixels, we calculate a slope of 0.03 pixels av- 
eraged over the geometrically corrected image (maximum 
slope of 0.99 pixels), while the average slope in the raw 
images is 2.54 pixels (maximum 3.43 pixels). Horizontal 
spectra are easier to extract and analyze using a vari- 
ety of tools. The orientation also allows simple filtering 
to remove cross-dispersion structure and isolate emission 
lines (e.g. Sec. 13. 2p . Finally, Apsis processing of the im- 
ages provides excellent cosmic-ray and hot-pixel removal 
and removes a small amplifier step (typically having an 
amplitude of a few electrons) often seen in WFC images. 

2.2. The data 

The HDFN field (RA = 12^36"^47.11% Dec = 
+62°13'11.9'') was observed by the ACS science team 
(program number 9301) for 2 orbits in the F775W (^775) 
filter, and 3 orbits each with the F850LP (2^850) filter and 
G800L grism as summarized in Tabled) Two exposures 
per orbit were obtained in order to facilitate cosmic-ray 
removal, and the telescope was dithered by 1 pixel in 
each axis between orbits. The individual CALACS pro- 
cessed G800L "FLT" frames were divided by the stan- 
dard F814W fiat-field image. Fine alignment of the in- 
dividual images was performed with Apsis. Apsis com- 
bined the individual exposures to make a single aligned 
image for each filter, F775W, F850LP, and G800L, using 
a spatial sampling of 0.05'' per pixel, as well as a detec- 
tion image that is the inverse sky variance weighted sum 
of the F775W and F850LP images. Because most ob- 
jects in grism images are rather elongated and faint. Ap- 
sis could not accurately register the G800L FLT images 
and determine the offsets. Instead it employed default 
shifts determined from the positions which are stored 
in the image headers. We used a modified version of 
Apsis to check the shifts in the G800L "CRJ" images 
from the STScI CALACS pipeline - these are pairs of im- 
ages combined to form a single cosmic-ray rejected image. 
There are three CRJ images for this dataset (one for each 
G800L orbit). Apsis matched 8-10 zero-order images of 
bright sources compared to the reference G800L CRJ im- 
age, yielding average shifts accurate to ~0.05 pixels in 
each axis. The resultant shifts matched the default shifts 
to 0.1 pixels. The final Apsis drizzle cycle was done to 
an output scale of 0.05" per pixel with interpolation per- 
formed using a Lanczos3 ke rnel. The Lanczos3 function, 
defined by Mefet^ (|2005f ). is a damped-sinc function. 
Application of it during drizzling results in better preser- 
vation of the noise characteristics and spatial resolution 
of the data t han the standa rd linear (square) interpola- 
tion kernel ( Mei et all [20051 ). The FWHM resolution of 
the final F775W and F850LP images was measured from 
direct measurement of stellar radial profiles and reported 
in Table [TJ For the G800L image we measure the spectral 



resolution from cross-dispersion cuts (five column sums) 
of the first order spectra of stars. These were fitted with 
a Gaussian profile, and the resolution taken as the av- 
erage FWHM of the fits. The resulting resolution of 2.1 
pixels corresponds to R = A/ A A 90 at A = 8500 A at 
the center of the field. 

We used th e fiux calibration curve given in 
IWalsh fc Pirzkall (|2004l ) to convert spectra to fiux units. 
As noted above, there is a ~ 10% variation in the fiux 
scale across the field when applied to data processed out- 
side of the standard aXe extractions from FLT frames. 
We employed a A calibration determined from WFC 
G800L images of Wolf-Rayet stars (which have strong 
bright emission lines) that were observed so as to fall on 
various positions on the WFC detector. The data and 
measur ement techniques employed are id entical to those 
used bv lPasquali. Pirzkal fc WalshI (|2003l ). but applied to 
the calibration data after drizzling them onto a rectified 
pixel grid with an output pixel scale of 0.05", the same 
as our data. The resultant A calibration is given as a 
quadratic polynomial as a function of column offset from 
the geometrically corrected direct image, with the poly- 
nomial's coefficients varying quadratically with position 
in the corrected direct images. 

3. FINDING EMISSION-LINE GALAXIES 

We have developed two methods for the semi- 
automated identification and classification of emission- 
line galaxies which we detail here. Here we define the 
term "emission-line galaxy" (ELG) to be a galaxy hav- 
ing line emission detected in our grism images. In prin- 
ciple, line emission in an ELG may be dispersed evenly 
throughout the galaxy. In practice it is usually confined 
to a small region, such as a nucleus or a knot. We use 
the term "emission-line source" ELS to denote a source of 
line emission that is distinct in position and wavelength. 
Effectively an ELS can be isolated as a distinct source 
in the grism image. Hence, an ELG with two knots each 
with only one detected line has two ELSs. However, if 
each of its knots has two distinct lines then there are four 
ELSs in the system. 

3.1. A: aXe selection 

Method A (for aX e) is very similar to that employed b y 
the GRAPES team (jPirzkal et al.l[2003 : IXu et al.lf2006h . 

The extractions are done using aXe with the calibra- 
tions discussed above (Sec. 12. 2p encoded into its config- 
uration file. The extractions are done from the Apsis 
processed G800L image which has a low order sky back- 
ground subtracted from it. No additional sky subtrac- 
tion was performed. The aXe inp ut catalog was derived 
from a SExtractor (|Bertin fc ArnoutsI Il996i) catalog of 
the detection image. Since we are starting with spectra 
of entire galaxies (and stars), any line emission we find 
can be attributed to the galaxy as a whole (making it an 
ELG) but not localized further. However, multiple ELSs 
can be discerned within an ELG if there is more than 
one emission line. We configured aXe to extract spectra 
down to the detection limit ~ 28.7 ABmag in the de- 
tection image. We go this faint to maximize the chance 
that we find faint "pure" ELGs - galaxies that have one 
emission line and no continuum. The grism observations 
have a similar exposure time and system throughput as 
the direct images, hence a pure ELG will have similar 
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count rates in the direct and grism images. Taking the 
differing exposure times and the rms noise level of the 
sky into account we calculate that the S/N ratio of a 
pure ELG would be similar in our detection and grism 
images. We caution that this condition may not hold for 
other datasets. 

The detection image catalog was processed to remove 
sources with m > 28.7 ABmag (too faint) and semi- 
minor axis size b < 0.8 pixels (too small, most likely 
image defects). aXe was configured to set the extraction 
aperture equal to 2.5 times the projected semi-major axis 
size, a of the objects. However, first we reset the size of 
sources having a < 2 pixels or 25.4 < m < 28.7 ABmag 
to a = 2, 6 = 2 pixels. Hence our smallest extraction 
aperture is 5 pixels wide. The final step before extract- 
ing the spectra was to mask the area within 8 pixels of 
the edges of the CCD chips in the grism image due to 
the number of false positives we found in preliminary 
runs with our code. 

One-dimensional fiux-calibrated spectra were ex- 
tracted with aXe. Automatic identification of "interest- 
ing" targets (ELG candidates) was performed by sub- 
tracting a smooth baseline spectrum and finding the 
sources with residuals having a peak S/N > 4. The base- 
line was constructed by median filtering the spectrum 
and then boxcar smoothing it, using a filter size of 19 
pixels in both steps. The spectrum of each candidate was 
displayed and classified as either (a) an ELG, in which 
case one or more Gaussian components were fitted to the 
peak(s) in the spectrum; (b) a "star", that is a source 
with strong broad absorption lines - our algorithm often 
mistakenly identifies the peaks between the absorption 
features as emission lines (the absorption lines sources 
are typically late type M or K stars alth ough we also 
found the two supernovae discussed by Blakeslee et al. 
[2003h . or (c) a spurious source. The sources classified in 
this manner are discussed in Sec. [H 

3.2. B: Blind grism image selection 

Method B (for blind selection) starts with the grism 
image and is designed to identify all detectable ELSs 
within galaxies. It is a "blind" selection in the sense 
that we do not require the a priori knowledge of source 
positions to find the ELSs. We find this to be very use- 
ful for two reasons. First, the ELSs we find are often 
confined to nuclei, off-center starbursts or strong H ll 
regions. Normal aXe extraction, as in method A, may 
dilute the line signal with "unnecessary" continuum flux 
or report the incorrect A for the line if it results from an 
off-center knot. This is because in aXe extractions, the 
A of each pixel depends on its offset from the major axis 
in the direct image; A errors may then occur for knots 
offset from the major axis. Since the flux scale depends 
on A, a flux error will also result. Second, method A can 
not flnd all possible pure emission-line sources. A pure 
ELG that emits at 5600A < A < 7100A wih be invisi- 
ble in our direct images since the the fllters we employ 
do not have signiflcant sensitivity at these wavelengths. 
However, the grism does, and hence we may still hope to 
flnd such sources, if they exist, in our grism images. 

Processing starts with high-pass flltering both the 
grism image and the detection image. This is accom- 
plished by smoothing with a 13 x 3 median fllter and then 
subtracting this smoothed image from the original. The 



long axis of the fllter is parallel to the image rows, that 
is, very nearly parallel to the dispersion direction. The 
flltering effectively removes most of the continuum in the 
grism image and much of the low frequency spatial struc- 
ture of the detection image, leaving compact ELSs in the 
grism image and galaxy nuclei, bars and knots in the 
direct images. The high-pass flltered grism images also 
contain the zero-order images, offset by ~ —115 columns 
from the direct images. These could be mistaken for 
ELSs. So, before searching for emission-line candidates 
we mask those that could contaminate our results. This 
is done by determining a linear coordinate transforma- 
tion^ between the detection image positions and that of 
the the zeroth-order images as well as a mean flux ratio. 
In terms of count rate, objects in the F775W (F850LP) 
image are on average 32 (21) times brighter than their 
zeroth-order counterparts in the grism image. We ap- 
ply the appropriate flux ratio to the detection image to 
locate pixels that would have a flux equal to or greater 
than the sky noise in their zeroth order. The coordinate 
transformation is used to determine their location in the 
grism image. This pixel distribution is grown by a radius 
of three pixels by convolving it with a circular top-hat 
function. The resultant masked pixels are set to 0.0 in 
the flltered grism image. The total usable area of the im- 
age is then 1659972 pixels or 11.83 arcmin^. Using this 
masking, about 60 spurious sources are excluded from 
the source catalogs (discussed below), while only 0.13% 
of the otherwise good area of the image is masked out. 
Hence, the masking is very efficient at removing spuri- 
ous sources yet unlikely to remove many real ELSs from 
the grism image. Zero-order images may survive near 
the image edges where the direct image falls outside the 
fleld of view of our detection image. This condition is 
easily tested. Stars and very compact sources also re- 
main in the high-pass flltered images because they are 
sharper than the smoothing box cross-dispersion width. 
However, they are easily recognized and flagged in the 
classiflcation stage. 

We use SExtractor to flnd sources in the masked and 
flltered grism image. By experimentation, we found that 
setting SExtractor parameters DETECT.THRESH and 
ANALYSIS_THRESH to 1.15 and DETECT_MINAREA 
to 3 was suflicient to flnd most obvious compact line emit- 
ters visible by eye without introducing large numbers 
of spurious detections. We removed from this catalog 
sources with output parameters ELONGATION (axial 
ratio) greater than 2.5 (likely spectral continuum resid- 
uals), BJMAGE less than 0.4 or FWHMJMAGE less 
than 1 pixel (likely residual cosmic rays or hot pixels), 
or FWHMJMAGE greater than 7 pixels (spurious since 
sources this large should have be missing from the high- 
pass flltered images). For each source we extract a region 
extending from —150 to +10 columns from its position 
in the grism image and having an extraction width Ay 
rows equal to 1.25 times its size projected onto the cross- 
dispersion axis: 

Ay = 1.25V(asin(9)2 + (6 cos (9)2. (1) 

Here a, b are the semi - major, minor axes AJMAGE, 
BJMAGE from SExtractor and is the position angle 
measured counter-clockwise from the -\-x axis deflned as 

9 Xout = a + bXin + cFin; Yout = d + eXin + /Fin 
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having constant row number in the pixel grid and di- 
rected towards increasing column number. The +x axis 
is close to, but not exactly the dispersion axis directed to- 
wards increasing wavelength. We set a minimum A?/ = 5 
rows for the extraction. This region is extracted from 
both the grism and direct filtered images, then the rows 
are summed to make ID-cuts. The region outside of a 
13 pixel box centered on the emission line is set to 0.0 in 
the grism cut to isolate the ELS. The grism and direct 
cuts are then cross-correlated to determine the x offset 
between the ELS and sources in the direct image. If the 
cross-correlation peak corresponds to the correct source 
in the direct image, then the x offset yields the source 
position in the direct images and hence a preliminary 
estimate of the wavelength of the line. 

Final measurements of the emission-line quantities are 
obtained from ID spectra of each knot extracted with 
aXe using the cross-correlation determined position, and 
an extraction aperture of 5 pixels. The emission-line 
properties are measured with Gaussian fits as in Sec. 13.11 
However, we use a peak S/N = 3 cut, lower than that 
employed by method A, since we find that method B can 
(usually) reliably find ELSs at this low of a significance 
level. Comparison of the fluxes of 19 single knot ELSs 
found by both methods show that lines are on average 
0.04 dex brighter (with a dispersion of 0.14 dex) when 
measured with method A compared to B. We consider 
this not to be a significant difference. 

Figure [T] shows an example of the image manipula- 
tion and detection process for this method. Since there 
are a number of ways this method can produce spurious 
results, this technique is applied in an interactive en- 
vironment. For each candidate ELS, the cutouts of the 
grism and direct images are examined, and used to assess 
whether it corresponds to a blemish in one of the images 
or a star (like method A, this technique is also adept 
at finding compact broad absorption-line sources). Line 
plots of the ID cuts and the cross-correlation spectrum 
are produced, as is the auto-correlation of the high-pass 
filtered grism spectrum with itself. The peaks in the 
cross-correlation are fitted with Gaussian profiles until 
the residuals have no peaks with S/N < 3. If there is 
more than one Gaussian component in the fit, the correct 
match is interactively selected using the 2D cutouts and 
line plots as a guide - an ELS typically corresponds to a 
high-surface brightness compact nucleus or knot with its 
2D line image resembling the high-pass filtered grism im- 
age in size and orientation. Spurious peaks in the cross- 
correlation typically can be identified (by eye) as having 
the wrong cross-dispersion position or not corresponding 
to a nucleus or knot. There remain cases, however, of 
more than one plausible direct counterpart to the ELS. 
This could be due to multiple knots in the direct image, 
or completely separate sources. We flag these ambigu- 
ous cases. The sources classified by this technique are 
discussed in Sec. IH 

3.3. Line identification and redshift 

Two emission lines are found in seven ELGs, allowing 
line identification and redshifts to be determined using 
the ratio of observed wavelengths, which is invariant with 
redshift. In three cases the two lines are rather close 
and clearly correspond to and the [O III] A4959,5007A 
doublet, which is blended at the grism's resolution. In 



two cases the ratio of wavelengths indicates that the 
lines are Ha and [Olll]. Note that one must be careful 
with this technique since Xiia/\oiu] = 1.3138 is close 
to Ah/3/A[oii] = 1.3041. A one pixel uncertainty in both 
line wavelengths could result an ambiguous line identifi- 
cation over the redshift range of interest. Our adopted 
identification in these two cases corresponds to previous 
spectroscopic and photometric z estimates (see below), 
indicating that our identifications are correct. In one 
case we identify the lines as RfS and H7 at z = 0.947, 
while in the last case we identify [Oil] and [NeIIl]3869A. 

The vast majority (39/46) of ELGs found contain only 
a single detected line. Identification of this line is a pri- 
mary, but difficult, task. The line may well be a blend 
at the low resolution of our data (e.g. Ho; and [Nil], the 
[Olll] doublet, and the [O Il]A3726,3728A doublet). Sim- 
ilarly, the resolution is not high enough to identify lines 
by profile shape (e.g. Lya). For these sources, our ap- 
proach is to use available redshifts as a first guess to each 
source's redshift and then determine which line identifi- 
cation agrees best with the guess. 

We use four sources for the redshift guesses, one spec- 
troscopic source and three photometric redshift sources. 
Multiple estimates are used to insure that all grism 
sources have at least one first guess redshift, and as 
a consistency check to determine how resilient the line 
identification is to the first guess source. In addition, 
the different redshift sources represent different choices 
of strategies, and investments in telescope time, labor 
and resources that may be employed to obtain red- 
shifts. We have four redshift sources which we now 
list with the number of matches to the 46 ELGs [in 
brackets]: (1) spectroscopic red shifts (zspec) from the 
compilation of the Hawaii group (|Cowie. et al .1120 04) [22 
match es] ; (2) photometric redshifts (^hot) from ICapalj 
(|2004l . hereafter C04) who use terrestrial U, B, V, 
i?, /, and HK photometry and derive 2:phot esti- 
mates from version 1.99 of the B PZ code described by 
Benitez (2000) [37 matches]; (3) iFernandez-Soto et al.l 
(1999, hereafter FLY99) who use -^450? ^06? 

and /814 photometry which th ey measure from the HDF 
data of I Williams et al.l (|l996l ) as well as g:rQund-ba sed 
J, and K photometry from iDickinsonl (|l998l ) to 
derive their %hot estimates [25 matches]; and (4) our 
own Zphot estimates which we derive from the on-line 
GOO DS ^485,1^06,^775,^850 ACS photomctry (release 
rl.lz. iGiavalisco et al.ll2004a[ ) [43 matches]. For the lat- 
ter we employ ed version 2. 0-alpha of the BPZ code de- 
scribed by Ben itez I (j2000l ) which reports results on up to 
three peaks in the z probability distribution. We assume 
that the emission lines seen are one of the following: Ha, 
the [Olll]4959,5007A doublet, the [O ll]3726,3729A dou- 
blet, or Lyo;, for which we adopt rest- wavelengths in vac- 
uum of Ao = 6564.6,4996.5,3728.7, 1216A, respectively. 
We adopted intrinsic flux ratios of ^5007/^4959 = 3.03 
for the [Olll] doublet (set by quantum mechanics), and 
^3728/^3726 = 1-3, Corresponding to an elect ron density 
Ue = 100 cm~^ (near the low density limit; lOsterbrocS 
[T989h for the [Oil] doublet. 

4. RESULTS 
4.1. Data presentation 
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Properties of the 53 ELSs in the 46 unique ELGs found 
by our techniques are hsted in Table [21 The first column 
gives the ELG identification number. The sources are 
ordered by ^grism, larger ELG numbers indicate higher 
redshifts. The second column gives a coordinate based 
name, mostly from the GOODSN rl.lz catalog. Remain- 
ing columns are explained in the table notes. Postage 
stamp cutouts of the ELGS taken from the detection im- 
age are given in Fig. [2l Two stamps are given for each 
source: one from the plain detection image and the other 
from the high-pass filtered detection image. The former 
is ideal for identifying the field while the latter is well 
suited to show knots and other sharp structures not al- 
ways seen in the plain detection image. Comments on 
individual ELGs are given in the appendix. 

We found 32 ELSs in 30 ELGs using method A - two 
of these ELGs have two emission lines detected with 
method A. Method B reaped more ELSs - 49 in 39 unique 
ELGs. Five of these have two lines detected with method 
B. Two ELGs have two identified emission knots, and 
one has four emission-line knots. The object centers for 
the extracted spectra and the method of detection is in- 
dicated in the stamps: squares mark sources identified 
with method A and circles mark sources from method B. 
The line identification corresponds to [O III] in 26 cases, 
[Oil] in 13 cases. Ha in eight cases, H/3 in four cases, 
H7 in one case, and [NeIIl]A3869A in one case. There 
were four candidate "pure" ELGs - apparent emission- 
line sources with no direct counterpart. However, careful 
examination of the grism image and noise maps produced 
in the drizzling process show that they are all probably 
imperfectly removed cosmic rays and hot pixels. 

4.2. Properties of grism- selected ELGs 

Here we examine some basic measured properties of the 
sample. These include line wavelength, flux, and equiva- 
lent width, as well as continuum magnitudes, colors, and 
luminosities. Our intent is to give an overview of the 
properties of our sample, examine the extent to which 
they are set by selection effects, and contrast them with 
other distant galaxy samples. It is not an analysis of 
the astrophysics of the ELGs, nor do we address issues 
of cosmology and the evolution of galaxy populations, 
since the emphasis of this paper is on how to find ELGs. 

Figure [3] compares the wavelength distribution of our 
sample with the grism sensitivity curve. The distribu- 
tions from method A and method B selection are shown 
separately. Both distributions and the curve peak at sim- 
ilar A and have a relatively long red tail. However, there 
is deficit of detections of lines at A < 7000A compared 
to the sensitivity curve. One-sided Kolmogorov-Smirnov 
tests yield probabilities of 7% and 3% that the observed 
A distribution follows the grism sensitivity curves for 
method A and B selection respectively. So while the 
sensitivity curve may play a significant role in determin- 
ing which ELGs are selected, other factors, including the 
volumes accessed by the individual lines, the luminos- 
ity function and large scale structure also effect the A 
distribution. 

We compare our line flux distributions in Figure [H 
The observed distributions are similar for the two selec- 
tion techniques, although there is a hint that method 
B is finding more faint lines than A. The peak in the 
distributions is at F^ne ^ 4 x 10~^^ergcm~^ s~^. For 



comparison we show the Fiine distribution from the STIS 
Parallel Survey f SPS iTeplitz d~^[2003i] lbl) . The SPS is 
similar to our survey in that it was obtained with HST 
using observations with a range of exposure times (typ- 
ically a few thousand seconds, i.e. similar to that used 
here), albeit with an instrument, the Space Telescope, 
Imaging Spectrograph, which has a much smaller collect- 
ing area and lower throughput than our observing con- 
figuration. We see that our observations typically reach 
about four times deeper in line flux, probably because of 
the improved sensitivity. The peak in our F^ne distribu- 
tions corresponds to a lOa detection within a 5 by 5 pixel 
pixel box for a line having A at the peak of the sensitivity 
curve. However, the A distribution is broad, and many 
lines are detected off the peak in the sensitivity curve 
(cf. Fig. [3]). The distribution of measured signal-to- noise 
ratio (S/N) is shown as insets in Fig. [H The S/N his- 
tograms show that we start to lose lines dit S/N < 6. As 
noted in sections 13.11 and 13.21 our initial source selection 
requires a peak S/N = 4 and 3, respectively in ID spec- 
tra. We are not finding weaker detections because here 
we are showing S/N measurements of integrated line flux 
within a 5 by 5 pixel box. One could dig further down 
in flux by convolving the ID spectra with a matched 
filter, and then lowering the peak S/N detection limit, 
although this would require weeding out more spurious 
sources. 

Figure [5] shows the ^775 magnitude distribution of the 
ELGs in comparison to two other sample s in this field 
- the Hawaii group Zspec compilation of 'Cowi e. et al.l 
(|2004D and the %hot catalog of C04. Only the sources 
that match sources in our detection image are included 
in the histograms to insure that the comparison cov- 
ers the same area on the sky. For the grism sample, 
the median ^775 is 23.9 ABmag and the 75th percentile 
is Z775 = 24.9 ABmag. These are significantly fainter 
than ^775 = 23.1,23.8 ABmag respectively at these per- 
centiles for the Hawaii 2;spec compilation. Through ex- 
pending considerable effort and time (particularly with 
the Keck telescopes) they were able to determine spec- 
troscopic redshifts for over four times as many sources 
(210) than the number of ELGs we find. Nevertheless, 
with a modest three orbit exposure we determined grism 
redshift estimates for 23 sources with no previous z^pec 
estimates. 

The faint apparent magnitudes correspond to low lu- 
minosities compared to typical galaxies. This is shown on 
in Fig. [6] which plots histograms of absolute magnitude in 
the rest-frame B band, Mb for our ELG sample, divided 
by line identification. We derived Mb by interpolating 
the broad band SED of the E LG host from the H ST 
GOODS optical photometry (G iavalisco et al.ll2QQ4af ) to 
determine the apparent ABmag at rest Aq = 4297A, the 
central A of the B filter. We used the apparent magni- 
tude in the nearest filter with data for cases where this 
A in the observed domain falls outside the range of cen- 
tral wavelengths of the ACS filters used by GOODS. The 
absolute magnitude was then calculated for our adopted 
Zgrism assuming Hq = 70 km s~^ Mpc~^, Q m = 0.3, 
and = 0.7 (following "CarrolLelalJ [l992|). Table E] 
presents the first quartile, median, and third quartile 
values of Mb and ^grism for our ELG sample, divided 
by emission line identification (ID). We use a thick gray 
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line to mark the absolute magnitude corresponding to 
m = 27 ABmag for an object at the redshift correspond- 
ing to the relevant line being found at A = 7500A (the 
peak throughput of the grism). This gives a crude in- 
dication of the limiting magnitude for each type of line 
emitter. It is not an absolute limit, since less luminous 
sources are still possible corresponding to finding some- 
what fainter apparent magnitudes (see Fig [5]) or bluer 
lines (hence closer sources). For comparison, we also 
mark the characteristic absolute magnitude at the knee 
of the luminosity function, M^, derived from various 
galaxy surveys in the literature as vertical broken lines in 
Fig. m The galaxy surveys used were chosen to sample 
similar redshifts as the ELG samples and to be broadly 
representative of the field galaxy population. For the 
Ha emitting ELGs we show derived from the Sloan 
Digital Sky Survey luminosity function (|Blanton. et al.l 
I2QQ3L = -20 31) and the Two Deg ree Field Galaxy 
Redshift Survey (|Madgwick et al.ll2QQ2l . M% = -20.56), 
while for the [Olll] and [Oil] emitters we use the K 
band selected rest frame B band GOODS luminosity 
function results at z = 0.46 (M% = -21.43) a nd 0.97 
(M^ = -21.45) respectively ( Dahlen et al.ll2005r ). 

In all cases the ELG samples have median Mb values 
below that of the field population at the same redshift. 
As one might expect, the difference between and the 
median Mb of the ELGs depends on the redshift and 
hence the lines identified in our survey. There is an in- 
teresting difference in the histograms at the bright end. 
There are no Ho; nor [Olll] emitters with Mb < —21, 
while seven of the eleven [O II] emitters are brighter than 
this value, with the most luminous, object #44, having 
Mb = —23.5. This difference may be due to sample size 
and relative volume: the volume available to [O II] emit- 
ters is 14 times that available to Ho; emitters and 2.7 
times that available to [O III] emitters. Hence if the sam- 
ples all had the same parent luminosity function, then 
the density of [Oil] emitters implies that we might ex- 
pect about three [Olll] emitters with Mb < —21 in our 
survey. A larger sample is needed to determine if the ap- 
parent deficit of high luminosity [O III] emitters is real. 
The surface density of the most luminous [O II] emit- 
ters is consistent with what we know about the luminos- 
ity function a t the redshifts we sample. For example, 
iDahlen et all (l2005l) find galaxies with Mb < -22 have 
a surface density of 0.37 arcmin"^ from a Ks band selec- 
tion in the GOODS-S field, while the three galaxies we 
find yield a surface density of 0.25 arcmin"^. 

The rest frame equivalent width EW is effectively a 
measurement of the ionizing fiux relative to the under- 
lying continuum. The EW was calculated from the aXe 
spectra using 



where fx is the continuum fiux density, measured from 
the spectra. Caution must be used in interpreting EW 
values, especially at the high end, due to background 
subtraction uncertainties. For Ha, the EW gives an in- 
dication of the present star formation rate compared to 
the past average. The interpretation is less clear for the 
[Olll] and [On] EW. [Olll] is considered a less reliable 
tracer of star formation because of its metallicity depen- 
dence, while the underlying continuum for [O II] emission 



is likely to be dominated by A stars perhaps produced in 
the same event causing the [O II] emission. 

Figure [71 shows the EW distribution of our sample 
compared to a variety of galaxy samples out to mod- 
erate redshifts. Equivalent width statistics of our sam- 
ple, split by line ID, are compiled in Table [31 as are 
the corresponding statistics of the comparison samples. 
Figure [7^ compares all the lines we measure to two lo- 
cal samples: ELGs found with prism spectroscopy by 
the KPNO International Spectroscopy Survey (KISS) red 
surveys (Salzer et al. 2001; Gronwall et al. 2004), and Hi 
selected ELGs imaged with narrow-band filters for the 
Survey of Ion ization in Neutral Gas Galaxies (SINGG 
Im eurer et all [2006). The grism-selected ELGs have sig- 
nificantly higher EWs. The differences between the sam- 
ples may largely be due to selection, measurement or 
instrumentation differences. While both the KISS and 
our grism survey use slit less spectroscopy, the KISS sur- 
vey has a higher dispersion (24 A pixel"-^), and employs 
a filter that limits the spectral range to cover 800A, thus 
limiting the sky background. This makes it easier to de- 
tect lower EW systems. Because of the higher spatial 
resolution of our data, the ACS spectra have extraction 
aperture widths of typically 0.25'', more than an order of 
magnitude smaller than the 4''-5'' resolution of the KISS 
survey. This means that the continuum in the ACS grism 
spectra are less diluted by non line-emitting portions of 
the host galaxy. The SINGG survey uses a totally differ- 
ent technique - narrow-band images to isolate Ha and R 
band images for continuum subtraction accurate to a few 
Angstroms, allowing even lower EW values to be mea- 
sured. Consequently SINGG includes many low surface 
brightness and low EW systems. 

Figure [7)3-d splits our sample by line ID and compares 
the EW histograms to those from the Canada-France 
Redshift Survey (C FRS [Hammer et al.[ [l997) and (SPS 
[Teplitz et aLll2Q03al lbl) for Ha emission (panel b), [Olll] 
emission (panel c), and [Oil] emission (panel d). The 
CFRS is an / band selected terrestrial spectroscopic sur- 
vey, while the SPS is a HST slitless spectroscopic sur- 
vey of effectively random high-latitude fields. Both com- 
parison samples have numerous detections of the three 
lines of interest and samples that extend out to z ^ 1.5. 
Our survey typically finds higher EW values for all lines 
than found in the CFRS. Again, this could in part be 
due to instrumentation differences; the CFRS spectra 
were obtained through 1.75'' wide slitlets (jLe Fevre et al.[ 
[1995[ ). typically covering a large fraction of the galaxy, 
and about an order of magnitude larger than the extrac- 
tion apertures we use. The SPS data is closer in nature to 
ours. The STIS slitless spectral resolution is significantly 
finer than our data - the two pixel resolution element cor- 
responds to ~ lOA; allowing STIS to detect lower EW 
features. Despite that, the SPS EW distributions are 
broadly similar to ours. This suggests that the relatively 
high EWs seen by SPS and ourselves compared to the 
CFRS may result from the smaller projected aperture 
sizes afforded by the space based observations. 

We next consider the broad band optical colors of 
ELGs in order to assess whether they can easily be se- 
lected by color. The ^435 — ^775 versus 2775 — zgso two 
color diagram of the ELGs identified in this study are 
compared to all HDFN galaxies in Fig. [HI The colors 
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were chosen because they span the broadest A coverage 
of the WFC filter set, and because they do a reasonable 
job at separating sources by redshift. Only sources with 
Zphot < 1-62 are included to limit the comparison sample 
to galaxies which would have [O II] , or lines to the red, at 
measured A < 9750A where they could be detected with 
the grism. While there are very few ELGs in the red tail 
of the field galaxy color distribution, the ELGs are not 
particularly blue. This can be seen in Table [4] which tab- 
ulates the median, first and third quartile colors of the 
two samples. The samples are subdivided by 2:phot 

(field 

galaxies) and line identification (ELGs) so that colors at 
similar redshifts are compared. In general, ELGs do not 
show as much of a "red tail" to their color distributions 
but are otherwise similar to the general field population. 
This is seen by their third quartile colors which are dis- 
tinctly bluer for the ELGs compared to the field galax- 
ies, while the first quartile and median colors are within 
~ 0.2 mag of the field galaxies. There are two excep- 
tions - the Ha emitters are somewhat redder than low- 
redshift field galaxies in ^435 — ^775 and the [O II] emit- 
ters are distinctly bluer than the field galaxies, especially 
in ^435 — 2775. We attribute the latter to the shift in Aq 
of the filters with redshift. At z = 0.15,0.53,1.13 the 
^435 — ^775 color samples rest frame colors m378 — ^674, 
772284 — ^507, ^204 — ^363"^^ " hcucc wc are measuring rest 
frame ultraviolet colors for the [Oil] emitters, while the 
^435 — ^775 colors for the other line emitters is still largely 
an optical color. Presumably the ELGs are mostly star- 
forming galaxies. This suggests that the rest-frame opti- 
cal colors of ELGs are fairly normal, but the rest frame 
UV colors are blue. This is reasonable, and can be seen 
in the template spectra used in photometric redshift es- 
timates. When normalized to unity at A = 4000A there 
is relatively little variation in SED shape for A > 4000A 
but strong color variations with UV color getting bluer 
as spectral type becomes later. 

A better separation of the line emitters can be obtained 
by including broad band filters that extend further to 
the blue than WFC's filter set. Figure [9] shows that the 
(U — V) versus (B — z) diagram is particularly well suited 
for emission-line identification. We examined a variety of 
optical two color diagrams and found this to be the best 
at discriminating between line emitters. If data in only 
three filters can be obtained, then the {U — V) versus 
(V — z) diagram also provides reasonable color separation 
between the emitters of the different lines. 

5. LINE IDENTIFICATION AND REDSHIFT 

Figure [10] compares the redshifts from the grism 
data with the spect roscopic redshift measurements from 
ICowie. et al.l (|2004l ). Symbol shape and color are used 
to indicate the line identification, while closed and open 
symbols mark measurements with methods A and B re- 
spectively. The dispersion in Zgrism — ^spec about the 
unity line is 0.016 for method A (15 measurements) and 
0.009 for method B (20 measurements) after applying 
an iterative 3a rejection (1 measurement was rejected in 
each case). We attribute the lower scatter from method 
B as refiecting its superior accuracy in measuring A in 
off-center star formation knots. 

the number in the subscript gives the approximate central A 
in nm. 



The redshift accuracy depends on how secure the line 
identification is. For the relatively bright sources plotted 
in Fig. [TOl the identification can be considered secure in 
the cases where more than one line is identified or where 
the line is identified with Ha or [O II] . These two lines are 
relatively isolated and usually there is no other plausible 
line within the redshift uncertainty of the first guess. In 
these cases, the dispersion in ^rism— ^spec about the unity 
line is 0.010 for method A (10 measurements) and 0.007 
for method B (13 measurements). The line identification 
is not secure if there is only one line and we identify 
it as [Olll]. The line could also be H/3; however, by 
default [Olll] is adopted as the identification under the 
assumption that typically the [O III] doublet is stronger 
than fcf lHammer et al.l ll99/). There are 5 (6) ELGs 
with spectroscopic z and single lines identified as [Olll] 
in our data and found with method A (B); the dispersion 
of their residuals about the unity line is 0.025 (0.009) 

Single-line identification depends critically on the first 
guess redshift. We find that in ~ 90% of the cases with 
both spectroscopic and photometric redshifts that 2;phot 
is sufficiently accurate to get the correct line identifica- 
tion. Specifically, when ^^spec is not employed as one of 
the redshift guesses 1 out of 16 of the method A line 
identifications changes, while 2 out of 20 of the method 
B identifications change. Of course, one also has to 
be careful with the sp ectroscopic redshifts; as noted by 
iLanzetta. et al.l (|l999f ) the misidentification of spectro- 
scopic sources previously had resulted in large discrepan- 
cies between z^pec and ^hot in the HDFN. This is likely to 
be the case for ELG #37, which is the outlier in Fig. [TOl 
It has 2:spec = 0.341 in the Hawaii catalog which would 
imply that we might see Ho; at A = 8802A or [Olll] at 
A = 6657A. No features are seen near either wavelength. 
Instead we see a strong line at A = 7467A. This source 
has 2;phot = 0-97 from C04 and 2;phot = 0-92 from our own 
BPZ results. Using either of these produces a first guess 
line identification as [Oil] at %rism = 1-00, well within 
the expected redshift errors of both %hot estimates. 

The reliability of our ^grism estimates is similar to 
that seen for photometric redshifts, while the accuracy is 
much better. This can be seen by co mparing; photomet- 
ric and spectroscopic redshifts (from lCowie. inil[200l 
for the sources in the field. The dispersion about the 
unity line in %hot — ^spec is 0.073, 0.107, 0.082 for ^^phot 
estimates from C04, FLY99, and our BPZ results respec- 
tively. Here we adopted a |2;phot — ^specl > 0.32 rejection 
criterion and only considered sources with ^spec 

< 1.5. 

The rejection criterion corresponds to three times the 
dispersion in Zphot — ^spec from FLY99 calculated with a 
3(7 clipping, while the 2:spec limit is adopted to correspond 
to the observed Zgrism range of our survey. The number 
of sources rejected/used in these calculations are 11/160, 
3/103, and 12/107, respectively. Hence - 5% to 10% of 
photo-z estimates are significantly discrepant compared 
to ^^spec- This is similar to the reliability of our line iden- 
tifications using 2;phot as the first guess redshifts. More 
importantly, the 2:grism results are more accurate than 
%hot (have a lower dispersion in z — Zspec) by a factor of 
5 to 12. 

The above comparisons, of course, require a spectro- 
scopic redshift. As illustrated in Fig. [5] these correspond 
to the brighter galaxies. The mean ^775 is 22.7 ABmag 
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for the 22 ELGs we selected which also have ^^gpec- The 
24 ELGs without Zgpec have a mean ^775 = 25.0. For 
them the main reliability issue is the photometric error 
bars. The errors are typically larger for faint galaxies, 
and in six of these galaxies they translate into redshift 
errors large enough to allow alternate bright line identi- 
fications. Choice of Zphot source can also be an issue. In 
eight cases the line identification changes depending on 
which first guess redshift is used. In one galaxy without 
Zgpec, no line identifications are allowed in the range of 
allowed %hot from our one estimate of %hot- There are 
a total of 14 cases that are ambiguous in one or more of 
these ways, 12 of those do not have a ^^gpec estimate. The 
Z775 ABmag distribution of the sources with ambiguous 
line identification, are marked in gray in the top panel 
of Fig. [5l illustrating that the ambiguous identifications 
correspond to faint sources (mean ^775 = 25.46 ABmag). 
These cases are identified in Table [2] with alternate line 
identifications noted in the appendix. While some un- 
certainty remains for these objects, we emphasize that 
by selecting the line identification closest to the favored 
^phot and prioritizing our Zphot sources, we increase the 
probability that we have picked the right line. 

While one could hope that additional priors might re- 
move the ambiguity of the line identifications, we have 
not found a satisfactory mea surement to use. For exam- 
ple, iDrozdovskv, et all (|2QQ5l ) decide on line ID, in part, 
by looking at the size of the host galaxies. However, size 
alone is not a great indicator of redshift. This is demon- 
strated in Fig. [11] which shows the size versus redshift 
relationship of the objects in our field. There is little if 
any angular size evolution with redshift. 

Correlations between the luminosity and line ratios 
have also been suggested to us for improving identifi- 
cations. For example, the mass- met allicity relationship 
(|Tremonti et al. 2004) is in the sense that that galax- 
ies with low mass, and hence low luminosity, have low 
metallicities, resulting in high excitations and thus typ- 
ically higher [OlIl]/H/3 and [O III] /[Oil] ratios: low lu- 
minosity galaxies are more likely to be [Olll] emitters 
while high luminosity galaxies are more likely to be [O II] 
emitters. However, on its own luminosity is unlikely to 
be be useful in constraining identifications for (at least) 
three reasons. (1) Using G800L tends to select high EW 
systems which are more likely to be high-excitation, low 
met allicity systems; this should induce a bias towards 
[Olll] emitters. The calibrating sample would need to 
have similar selection effects as the grism ELGs. (2) The 
mass metallicity relationsh ip is known to evolve with red- 
shift (Savaglio et al.|[2"005f ) ; higher redshift galaxies of the 
same stellar mass have lower metallicities (again favoring 
[Olll] emitters). (3) The luminosity-redshift relationship 
goes in the wrong direction to remove the degeneracy. 
For example ELG #45, one of the cases with ambiguous 
line identification, has ^grism = 1.422 for our adopted 
[On] identification, yielding Mb = —21.3, brighter than 
the local which seems consistent with the [O II] iden- 
tification. If we adopt the alternative [O III] identification 
then Zgrism = 0.807 and Mb = —19.2, much fainter than 
M^, which seems to be consistent with an [Olll] iden- 
tification. Either combination seems plausible and the 
degeneracy is not broken. 

One may also consider using EW as a prior to decide 
between possible line identifications, particularly in con- 



junction with luminosity. In the local universe it is very 
rare to have EW{[Ou]) > 100 A; less t han 2% of the 
prism selected Ha ELGs in the sample of 'Galle go et al.l 
(|l996l ) meet this condition. Such EW values are more 
common from [O III] emission. When seen in [O II] , the 
source typically has a low (fainter than M^) luminosity; 
Mb ^. —20 (when converted to our adopted co smology 
iPerez-Gonzalez et al.|[2QQQt iGallego et al.lll996[ ). Using 
EW as a prior would cast suspicion on the seven single 
line [0 11] galaxies in our sample with such high EW val- 
ues. Such scrutiny is warranted since six of these seven 
have some ambiguity in the line identification (as noted 
in the appendix); in four of those cases an [Olll] iden- 
tification is allowed depending on which 2;phot is used 
as a first guess redshift. If we were to use luminosity 
as well in the prior, changing the [O II] identification to 
[Olll] in the ambiguous cases if £;V^([Oll]) > lOOA and 
Mb < —20.0, then two sources (#37 and #45) would be 
effected. 

While the identifications of the high EW [Oil] emit- 
ters deserve some skepticism, at this time it would be 
inappropriate to apply an EW prior even in conjunction 
with luminosity. There are three reasons for this assess- 
ment. (1) As noted in Sec. l4.2l the continuum levels used 
to determine EW are prone to large background subtrac- 
tion uncertainties hence the accuracy of high EW mea- 
surements are typically low. In this exploratory study 
we have not calculated the uncertainties in the contin- 
uum level which would need to be done to properly ap- 
ply a prior. (2) High £^W([Oll]) v alues have a prece- 
dence in the more distant universe: iHogg et all (fT998h 
find two [On] emitters with E W > lOOA, while the re 
are 22 such sources in the SPS (jTeplitz et al.l l2003al lbl) 
(these studies also note the problem in continuum deter- 
mination). All these cases have z^pec ^ 0.5. We cau- 
tion that the majority of the cases found by the SPS are 
also single line sources, and thus one might also be sus- 
picious of their proper identification. However, two of 
their high EW [Oil] emitters have additional lines that 
secure their identification. (3) Finally, we note that there 
is strong evidence for evolution in EW, with ^W^([Oll]) 
increasing with redshift (es pecially for z ^ 0.9), e ven for 
the most luminous galaxies [Hammer et all (jl997[ ). Rea- 
sons (2) and (3) indicate that an EW prior based on 
the local universe may not be appropriate at the red- 
shifts we are dealing with. We conclude that the high 
£^l^([Oll]) emitters need more scrutiny to confirm their 
reality. This should include a more careful determina- 
tion of £^VK([Oll]) using an improved error analysis, as 
well as follow-up spectroscopy to detect additional lines 
and confirm the line identification. If the high incidence 
rate of sources having £^VK([Oll]) > lOOA is confirmed, 
it would be further proof of strong redshift evolution in 
the star forming properties of galaxies. 

6. SUMMARY AND DISCUSSION 

We have shown that a modest expenditure, three or- 
bits, of ACS WFC grism time with the HST pointed at 
a "blank" high latitude field results in the detection of 
dozens of ELGs out to z ^ 1.5. Here we found 46 ELGs 
in the HDFN yielding a surface density of 3.9 ELGs per 
square arcmin. The "blind" grism selection technique 
(method B) results in significantly more sources and bet- 
ter redshift accuracy. We attribute this to its ability 
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to isolate individual emission-line knots within galax- 
ies. The aXe selection technique (method A) relies on 
an initial catalog of sources, and hence is effectively a 
broad band (i.e. usually continuum) selection technique. 
While it often misses objects where emission arises from 
a knot, it is adept at picking up line emission confined 
to a compact nucleus, which the "blind" technique can 
miss. Hence, the two techniques are complimentary. 

The ELGs found are most frequently [Olll] (or H^^) 
emitters at z ~ 0.4 to 0.9. Ho; and [Oil] emitters are 
also found, but are less common because of the smaller 
volume for the former, and the limited depth of the ob- 
servations for the latter. The ELGs represent a small 
fraction of the field population. There are 647 galaxies 
within the field of our observations having Zphot ^ 1-5 in 
the Capak (2004) catal og, while 186 galaxies have spec- 
troscopic z < 1.5 in the lCowie. et al.l (|2004l ) compilation. 
While grism selection of ELGs does not result in a sam- 
ple of the field that is in any way complete, the galax- 
ies selected do have interesting properties. In particular 
they tend to be low-luminosity high-EW systems. This 
suggests that they are experiencing an intense burst of 
star formation, or may contain an AGN. A high [Olll] 
EW suggests high excitation and low metallicity. Grism 
selection of ELGs may be a good means to locate the 
barely evolved building blocks of larger galaxies. 

Our results are consistent with deeper G 800L observa- 
tions reported b y the GRAPES team (Pirzkal et al.ll2004l : 
IXu et all I2006D who obtained 92 ks of G800L observa- 
tions of the HUDF (13 times longer than our HDF ob- 
servations) split into five epochs. They found 113 ELGs 
in a field having similar area, using an algorithm equiva- 
lent to our method A (although differing in some details); 
51 of these are brighter than our empirical line fiux limit 
^line = 3.0x 10~^^ergcm~^ s~^. This compares well with 
the 46 ELGs we find in the HDF. They also find [O III] 
sources to be the most frequently detected line, while the 
maximum F^ne is 2.2 x 10~^^ ergcm"^ s~^ in the seven 
Lyo; emitters they find, consistent with our non-detection 
of these sources. 

Optimal use of grism data to discover ELGs requires 
additional data. This is because the grism spectra typi- 
cally show only one emission line per object, hence iden- 
tification of the line from the low-resolution spectra is 
difficult, at least for the relatively short exposures used 
here. With longer exposures, often both the [O III] dou- 
blet and Iif3 lines can be seen in ELGs having z ~ 0.4 
to 0.9, hence the problem then becomes distinguishing 
between Ha and [O II] emitters. 

The accuracy of the line identification can be improved 
if there is a good "first guess" redshift for each source, 
either a spectroscopic redshift ^spec or a photometric red- 
shift ^phot- While the former produces the most accu- 
rate line identifications, there is typically little need for 
a grism spectrum of sources that already have a ground- 
based spectrum of sufficient S/N to determine a red- 
shift. Use of Zphot as the first guess requires additional 
photometric data from HST or other sources to derive 
the redshift. Without these additional data, or follow-up 
spectroscopy, it may be impossible to identify the line 
and hence determine the redshift, which seriously dimin- 
ishes the utility of the ELG discoveries. With a good 
Zphot first guess, lines can be identified with ~ 90% re- 
liability, similar to the 2;phot reliability, but resulting in 



redshifts accurate to ~ 0.01 (3000 kms~^). This is about 
an order of magnitude better than Zphot estimates and is 
sufficient for separating ELG members of rich clusters 
from the field. 

The requirement of additional photometry to obtain 
good emission-line redshifts amounts to a substantial ad- 
ditional investment of time and labor. In Sec 14.21 we ex- 
amined color-color diagrams that are useful for sorting 
the ELGs by line identification. In general filter combi- 
nations that span the full optical range seem to be the 
best suited for this purpose. The U — V versus B — z 
diagram provides the best discrimination, but requires 
wide field U band data which is hard to obtain from the 
ground, and impossible to obtain with ACS + WFC. The 
^435 — ^775 versus Z775 — 2^850 diagram does a reasonable 
job at separating the [O II] emitters from the [O III] and 
Hq; sources. Most of the discrimination comes from the 
^435 — ^775 color which is the single color best suited for 
line discrimination from the ACS + WFC filter set. How- 
ever, as shown in Fig. [8] it does not discriminate well be- 
tween Hq; and [O III] emitters. For that one needs to have 
a filter as far as possible towards short wavelengths so as 
to sample the rest frame UV at modest redshifts. An 
efficient solution of this issue at HST resolution would 
require the installation of WFC3. Until then, terrestrial 
{U — V) gives the best single color discriminator between 
line identifications. 

The reliability of the line identification decreases with 
decreasing brightness, as the increasing photometric er- 
rors can result in ambiguous line identification. The in- 
vestment of direct imaging time required to beat down 
the photometric redshift errors is larger than the time 
spent on the grism imaging; we spent three HST orbits 
imaging with G800L, while the GOODS direct images 
(^435 - 3 orbits, Vqqq - 2.5 orbits, 2775 - 2.5 orbits, zgso - 5 
orbits, at each pointing), as well as the original HDF im- 
ages and ground-based imaging were used to determine 
the first guess photometric redshifts. Even then, more 
than half of the ELGS we detected with ^775 > 24.5 AB- 
mag have some uncertainty in their line identification. Of 
particular concern are the faint [O II] identifications that 
have apparently large EW > lOOA. These could signify 
strong redshift evolution in the st ar forming population , 
as indicated by other studies (e.g. iHammer et aDll997l ). 
or could be (in part) spurious due to contamination of 
misidentified [O III] emitters, or large continuum place- 
ment errors. The HDFN is one of the best studied deep 
field, yet we still face these issues because of the faintness 
of the ELGs. 

We conclude that ACS G800L grism data with minimal 
direct images provides a useful means of locating ELGs, 
but without additional data, can provide only a limited 
interpretation of the nature of the sources. Photometric 
redshifts from broad band imaging can improve the re- 
liability of line identifications. However, these data are 
also expensive to obtain, and the results are still likely to 
be ambiguous for the faintest ELGs we can detect. Se- 
cure redshifts for these still require ground-based spec- 
troscopy. Fortunately the line fiuxes are easily within 
the reach of the current generation of 8m class telescopes. 
For example the Gemini Mult i- Object Spectrographs can 
detect an emission line with a fiux of 1.5 x 10~^^ erg 
cm~^ s~^ (one tenth of our limiting line fiux) and EW of 



Emission-line sources in ACS grism data 



11 



lOA at S/N ~ 5, in a single 900s exposure. An efficient 
strategy for finding and characterizing ELGs would then 
be to observe with ACS and the G800L grism to find 
the emission-line sources, employing broad band images 
(in say F814W and F606W) of a similar depth to locate 
the corresponding galaxy; then following up with ground 
based spectroscopy to secure the redshift and identify 
additional lines. 
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APPENDIX 
COMMENTS ON INDIVIDUAL SOURCES 

Here we present notes on ELGs with multiple ELSs (lines and/or knots), those with striking morphologies, and cases 
where the line identification is in some way ambiguous. For each source we list the ELG id number, the corresponding 
long name, and grism-redshift (in parenthesis). 

#1 - GOODSN J123641. 63+621132.1 (0.098): The lowest redshift galaxy shows Ra emission from a bright 
knot, pinpointed with method B, 1.07^' to E of the spiral galaxy core. This knot combined with a fainter knot to the 
east of the core results in a weak Ho; detection with method A. 

#2 - GOODSN J123644.75+621157.4 (0.124): The adopted identification of the single line in this source as 
Ha yields a redshift closest to the C04 Zphot- However, the error bars in C04 and our Zphot also allow an identification 
of [O III] , while the 

%hot from FLY99 only allows [O III] as the identification. 

#3 - GOODSN J123633. 16+621344.0 (0.126): The adopted Ra line identification yields ^gnsm closest to Zphot 
from C04. However, the error bars from C04 also allow line identifications of [Olll], or [Oil], while our BPZ results 
allow Ha and [O III] as the line identification. 

#5 - GOODSN J123648. 30+621426. 9 (0.136): Lopsided spiral, with a bright compact nucleus and bar. The 
spectrum extracted with method A shows weak Ha on top of a strong continuum. With method B, four Ha emitting 
knots are identified with lines at a consistent A (average A = 7456A). A fifth ELS appears to have emission at a 
discrepant A perhaps due to misidentification of the emitting knot in the direct image. 

#7 - GOODSN J123658. 06+621300. 8 (0.308): Inclined disk galaxy with a prominent knot offset by 0.35''. 
[O III] and 11(3 emission arise from the galaxy nucleus. 

#8 - GOODSN J123626.57+621321.2 (0.318): Our BPZ analysis provides the only Zphot source. The Zphot 
error bars allow the single line to be identified as Ha, [Olll], or [Oil]. The adopted [Olll] identification is closest to 
the third strongest 2;phot probability peak from BPZ. 

#9 - GOODSN J123650. 82+621256. 3 (0.319): Edge-on disk galaxy with two emission-line knots separated by 
0.96'' bracketing the nucleus. With method B, [Olll] is clearly visible in both knots. Ha is clearly detected in the 
southern knot, but is just below the detection limit in the northern knot. No emission lines are detected with method 
A, because the extracted spectrum does not fully contain the knots. 

#10 - GOODSN J123646.59+621157.5 (0.341): The error bars on the C04 Zphot allow [Olll] or [Oil] as the 
line identification for the single detected line. The closest ^grism niatch is with [O III] which is also consistent with the 
FLY99 and BPZ 2:phot analysis. 

#14 - GOODSN J123637.56+621240.4 (0.445): This galaxy is apparently interacting with GOODSN 
J123637.64+621241.3, which we also detect as an ELG(see below). We find one line, [Olll], with method B. 

#15 - GOODSN J123637.64+621241.3 (0.446): The dominant system in the pair with #14. Method B detects 
two lines identified as [Olll] and Ha, while only one low EW line, [Olll], is identified with method A. 

#17 - GOODSN J123657.30+621300.0 (0.465): A barred spiral with arms forming a pseudo ring. Line emission 
originates in an arm H ll region about V from the galaxy center. 

#20 - ACS J123636. 58+621336. 8 (0.478): This source is not present in the GOODSN rl.lz catalog, probably 
due to its faintness. Our own measurements of this source from data combining the GOODS images and other ACS 
images of the field yield photometry [^435, V^eoGjJiz^, ^85o] = [28.16 ± 0.14, 27.69 ± 0.08, 26.77 ± 0.06, 27.72 ± 0.15] 
ABmag, while measurements from the Capak (12004) images yield [U, B, V, R, I] = [29.16 ± 0.55, 28.22 ± 0.30, 
28.00 ± 0.28, 28.11 ± 0.34, 27.7 ± 0.6] ABmag and z > 26 ABmag (a non-detection). The relative brightness in ^775 is 
likely due to the single bright line we observe at A = 7386A whose flux (Table [2]) is consistent with dominating ^775. 
Using this photometry, BPZ yields a best %hot = 0-62 consistent with the line being [Olll]. However a second peak 
in the probability distribution at 2:phot = 0.1 means that an identification as Ha can not be ruled out. Because this 
source was not in the GOODSN rl.lz catalog, it was excluded from the statistics given in Table |4] (below). 

#22 - GOODSN J123655. 58+621400. 3 (0.551): Emission corresponds to a knot above the plane of an edge-on 
disk galaxy. 

#24 - ACS J123657.48+621212.0 (0.555): At first blush, this object appears to be the nucleus of a dwarf galaxy 
being shredded by an interaction with its neighbor GOODSN J123657.49+621211.2 projected 1.46'' towards SSE. 
However, that source has 2;spec = 0.669 (C04). Our redshift is from a single line detected with both methods identified 
as [Olll]. Even if the line were H^^ at 2: = 0.600, or H7 at z = 0.790, the redshift would be significantly discrepant with 
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its neighbor. Hence, the apparent interaction may be spurious and the sources an unrelated chance projection. 

#26 - GOODSN J123654.39+621434.7 (0.573): The nucleus of this modestly inclined disk galaxy shows [Olll] 
and 11(3. The emission lines are found with method A, not B, probably because the line emission is centered on the 
compact nucleus. 

#27 - GOODSN J123645. 53+621330.2 (0.670): The large error bars on the C04 Zphot allow the single line to 
be identified as Ha, [Olll] or [Oil]; the closest %rism match is the adopted [Olll] identification. This is consistent with 
the smaller error bars on the 2;phot from BPZ, while none of the brightest likely lines match the ^^phot from FLY99. 

#28 - GOODSN J123636.47+621419.1 (0.684): The off-center knot in this smah galaxy shows a broad, bright 
emission line well fit as [O III] blended with H/3. 

#30 - GOODSN J123647.24+621134.7 (0.717): BPZ is the only source of Zphot for this faint source, and yields 
allowed redshifts in the range 1 < %hot ^ 2. However, none of the typical bright lines can match this range and the 
observed A. The adopted [O II] identification corresponds to the line closest to the allowed range. 

#33 - GOODSN J123644.17+621430.5 (0.858): The large C04 Zphot error bars for this source allows the single 
line to be identified as [Oil] as well as the preferred [Olll]. The 2;phot from BPZ is consistent with our adopted [Olll] 
line identification. 

#34 - GOODSN J123652.97+621257.1 (0.943): The preferred Zphot from FLY99 identifies the single line as 
[On], while the BPZ 2:phot indicates the line may be [Olll]. 

#35 - GOODSN J123636. 63+621347.1 (0.947) This bright compact galaxy is detected in both H^ and H7. 
It is the only H7 source in the sample . The fiux ratio Fnp/Fn^ ^ 1.9 is close to the expected case B ratio of 2.1 
(for Tie = 100 cm~^ and Tg = 10^ K Dopita & Sutherland 2003). The lines are detected with method A but not 
with method B, probably because the line emission is centered on the compact nucleus. The case for interaction with 
GOODSN J123636.85+621346.2, a larger but slightly fainter lopsided spiral 1.81'' to ESE, is strong since one of its 
spiral arm seems to be connected to GOODSN J123636. 63+621347.1 in the high-pass filtered direct image, reminiscent 
of the M51/NGC5194 system (see Fig [2]). However, the spectroscopic redshift of the spiral is significantly lower (0.846, 
C04) casting some on this inference. The ^phot 

from C04 are consistent with these line identifications, while our BPZ 

%hot is too low. 

#36 - GOODSN J123649. 35+621155.4 (0.954): Compact galaxy with two lines arising in the nucleus identified 
as [On] and [Ne III] A3869A. The measured line ratio i^[Are///] /^[o//] ^ 0-38 indicates a high excitation: it corresponds 
to the 80th percentile in this ratio for the local ELGs which display both lines in the catalog of Terlevi ch et al.l (jl99lf ). 

#37 - GOODSN J123649.47+621456.9 (1.003): This sourc e has the lar g est d iscrepancy between z^pec and 
Zgrism in Fig.[TOl The original source for the reported z^pec = 0.341 is lCohen et al.l (120001). wher e the spectrum is given 
a quality code of 5: "one emission line only, reality uncertain, assume 3727" (|Cohen et al.lll999l ). that is, a single weak 
line in the spectrum. Both our BPZ results and those of C04 indicate very similar Zphot* 0.921, and 0.970, respectively. 
While there is no good match at z ~ 0.34 for the line we detect at A = 7469A, adopting either Zphot estimate for our 
first guess redshift identifies the line as [Oil] at our adopted ^grism = 1.003. Using this redshift, then the line found 
bylCoheneEaD (|2000D may be [Oil] 2470 A or [NeIv]2423A, if the weak line they found is real. 

#40 - GOODSN J123645.46+621357.3 (1.073): The preferred ^phot from FLY99 yields the single line identifi- 
cation of [On], while the BPZ 

%hot indicates the line may be [O III] . 

#42 - GOODSN J123653.51+621141.4 (1.263): The C04 Zphot yields our adopted [Oil] line identification. The 
%hot from BPZ is lower and outside the C04 error-bars, but does not allow any of the standard bright line guesses to 
correspond with the observed line A. 

#44 - GOODSN J123652. 77+621354.7 (1.346) The brightest knot of this chain galaxy has one line identified 
as [On]. 

#45 - GOODSN J123642. 55+621150. 3 (1.422): The adopted [O II] identification for the single line of this source 
is consistent with 2;phot from C04 and FLY99, while BPZ prefers a lower 2;phot consistent with an [O III] identification. 

#46 - GOODSN J123648.48+621120.7 (1.424): BPZ is the only Zphot source for this faint source; the error 
bars allow an [O III] identification as well as the adopted [O II] . 
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TABLE 2 
Emission line sources 



ELG # 




Name 


^775 




-^phot 


S 


•2'grism 


ID 


EW 


log(Fiine) 


M 


Notes 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(11) 


I 


vJ Vv Vv kJ 1 > 




19.78 


0.089 


0.090 


C 


0.098 


Hq; 


127 


— 15.52 


BA 


b 


2 


VJ Vv kJ 1 > 


T1 9Sfi44 7^-KS91 1 "^7 4 


25.06 




0.210 


c 


0.124 


Ha 


223 


— 16.43 


B 


cd 


3 


VJT W W kJ 1 > 


T1 2SfiSS 1 fi-kfi21 S44 D 


25.31 




0.220 


c 


0.126 


Ha 


155 


— 16.59 


BA 


d 


4 


VJ W J_v kJ 1 > 


T19Sfi4fi '^'^-KS914n7 Q 


23.95 


0.130 


0.160 


c 


0.128 


Ha 


115 


— 16.33 


BA 




5 


VJJ V_y V_y J_y kJ 1 > 




19.01 


0.139 


0.130 


c 


0.136 


Ha 


62 


— 15.52 


BA 


b 


6 


OOODSN 


J123651. 72+621220. 5 


21.62 


0.300 


0.320 


c 


0.302 


¥La 


58 


— 16.20 


BA 




7 


OOODSN 


T1 2Sfi58 06+621 SOO 8 


22.45 


0.319 


0.310 


c 


0.308 




9 


— 16.68 


B 


a 


















OIII 


8 


— 16.78 


B 




8 


VJ -L/ kJ i > 


Tl 2S626 57+621 S21 2 


26.56 




0.483 


B 


0.318 


OIII 


6418 


— 16.68 


A 


d 


9 


vjr J-V »j 1 >i 


712^650 82-1-621256 3 


22.61 


0.319 


0.310 


C 


0.319 


OIII 


250 


— 16.00 


B 


ab 


















Ha 


169 


— 16.30 


B 




10 


vJ J-V »J 1 N 


71 23646 5Q-I-621 1 57 5 


25.77 




0.440 


C 


0.341 


OIII 


248 


— 16.44 


BA 


(J 


11 


OOODSN 

VJ V>' V^' J--' kJ 1 N 


7123653 10-1-621438 4 


24.53 




0.430 


C 


0.376 


OIII 


123 


— 16.49 


B 




12 


vjr vy J-V kj 1 N 


71 23628 76-1-621 335 8 


25.99 




0.559 


B 


0.427 


OIII 


204 


— 16.56 


BA 




13 


ACS 


7123632 6Q-I-62123Q 1 


23.15 


0.458 


0.420 


C 


0.433 


OIII 


65 


— 16.46 


B 




14 


OOODSN 


J123637.56+621240.4 


22.23 


0.457 


0.153 


B 


0.445 


OIII 


12 


— 16.58 


B 




15 


OOODSN 


7123637 644-621241 3 


21.08 




0.436 


B 


0.446 


OIII 


81 


— 16.09 


BA 


a 


















Ho; 


155 


— 15.99 


B 




16 


OOODSN 


7123650 7Q+621221 7 


24.71 




0.440 


C 


0.450 


OIII 


254 


— 16.27 


BA 




17 


OOODSN 


7123657 30+621300 


21.40 


0.473 


0.420 


C 


0.465 


OIII 


89 


— 16.41 


B 




18 


OOODSN 


7123644 10+621248 2 

U X Zj <X VJtit: . X ?_/ | Zj X ZjtiO . Z^ 


21.62 


0.555 


0.540 


C 


0.476 


OIII 


20 


— 16.05 


A 




19 


OOODSN 


7123637 76+621235 6 


23.86 




0.480 


C 


0.477 


OIII 


172 


— 16.36 


B 




20 


ACS 


7123636 58+621336 8 


26.77 




0.620 


C 


0.478 


OIII 


315 


— 16.45 


BA 


Qg 


21 


OOODSN 


J123645. 24+621 108. 9 


23.41 


0.513 


0.580 


C 


0.505 


OIII 


162 


— 16.14 


BA 




22 


OOODSN 


7123655 58+621400 3 


24.17 


0.559 


0.590 


C 


0.551 


OIII 


147 


— 16.44 


BA 




23 


OOODSN 


J123645. 32+621 143. 2 


23.85 


0.557 


0.570 


C 


0.552 


OIII 


198 


— 16.15 


BA 




24 


ACS 


71 23657 48+621 21 2 

tJ X ZiOvJO 1 • t:0^^ vJZi X Zi X Zi . U 


23.29 




0.720 


p 


0.555 


OIII 


26 


— 16.59 


BA 




25 


OOODSN 


7123644 75+621144 1 


24.87 




0.670 


C 


0.562 


OIII 


91 


— 16.54 


B 




26 


OOODSN 


7123654 3Q+621434 7 


22.26 


0.577 


0.690 


C 


0.573 


OIII 


13 


— 16.50 


A 






















11 


— 16.54 


A 


a 


27 


OOODSN 


7123645 53+621330 2 

O -L ^kJKJ^kJ • 1 vJZ;X?_lOW»Z; 


25.38 




0.540 


c 


0.670 [0177] 243 


— 16.47 


BA 


cd 


28 


GOODSN 


J123636. 47+621419.1 


24.30 




0.700 


c 


0.684 0777 341 


— 15.81 


BA 


a 


















H/3 


115 


— 16.47 


B 




29 


OOODSN 

vJ J-V k_J 1 >i 


71 23646 Q6+621 1 33 


24.27 




0.700 


c 


0.685 


OIII] 198 


— 16.12 


BA 




30 


OOODSN 

VJT vv vy kji > 


71 23647 24+621 1 34 7 


27.10 




1.518 


B 


0.717 


Oil] 


121 


— 16.24 


A 


f 


31 


OOODSN 

V-T V-/ V-/ -I— / k_/ 1 > 


J123629. 72+621329. 9 


22.83 


0.746 


0.700 


C 


0.737 


OIII 


28 


— 16.15 


A 




32 


OOODSN 

VJT Vy W -L/ 1 > 


7123642 20+621420 Q 

U X Zj <X vJt: Zj . Zj ?_/ | U Zj X t: Z; ?_/ . ?_/ 


23.77 




0.850 


C 


0.841 


OIII 


287 


— 16.04 


B 




33 


OOODSN 

VJT W Vv -L/ 1 > 


7123644 17+621430 5 


24.66 


0.863 


0.925 


B 


0.858 


OIII 


175 


— 16.23 


BA 


cd 


34 


OOODSN 


7123652 Q7+621257 1 


25.63 




0.800 


F 


0.943 


[Oil] 


182 


— 16.32 


BA 




35 


OOODSN 

VJT vv vy -1—-' kj J- ^ 


7123636 63+621347 1 


21.44 


0.962 


0.386 


B 


0.947 


H7 


7 


— 16.34 


A 


a 


















R(3 


14 


— 16.07 


A 


a 


36 


GOODSN 


J123649.35+621155.4 


23.41 


0.961 


1.115 


B 


0.954 


[Oil] 


32 


-16.39 


BA 


a 
















[Nel77] 13 


-16.81 


B 


a 


37 


GOODSN 


J123649.47+621456.9 


24.11 


0.341 


0.970 


C 


1.003 


Oil 


107 


-16.69 


B 


ce 


38 


GOODSN 


J123654.45+621152.8 


24.40 




1.040 


C 


1.017 


Oil 


36 


-16.63 


B 




39 


GOODSN 


J123658. 30+621214.5 


23.33 


1.020 


0.970 


C 


1.026 


Oil 


40 


-16.32 


B 




40 


GOODSN 


J123645.46+621357.3 


26.17 




0.920 


F 


1.073 


Oil 


425 


-16.37 


BA 


c 


41 


GOODSN 


J123643.42+621151.9 


23.12 


1.241 


1.200 


C 


1.237 


Oil 


22 


-16.55 


BA 




42 


GOODSN 


J123653. 51+621141.4 


23.89 




1.490 


C 


1.263 


Oil 


46 


-16.35 


BA 


c 


43 


GOODSN 


J123644.98+621240.0 


24.07 




1.170 


C 


1.337 


Oil 


578 


-16.20 


A 




44 


GOODSN 


J123652.77+621354.7 


22.73 


1.355 


1.440 


F 


1.346 


oil 


41 


-16.20 


B 




45 


GOODSN 


J123642.55+621150.3 


24.97 




1.600 


C 


1.422 


Oil 


228 


-16.57 


B 


c 


46 


GOODSN 


J123648.48+621120.7 


26.87 




1.141 


B 


1.424 


oil 


198 


-16.49 


B 


d 



Note. — Column descriptions: (1) ELG catalog number (this work). (2) Names preceded by GOODSN are the lAU specified name 
from the GOODS release rl.lz (Giavalisco 2004a,b). Names preceded by ACS could not be matched with the GOODSN catalog. (3) 2775 
are SExtractor magnitudes of our images through the F775W filter in the ABmag system. (4) The spectroscopic redshift 2;spec taken from 
the Hawaii group compilation (Cowie et al. 2004). (5) and (6) are the adopted photometric redshift 2;phot its source: C - Capak (2004); 
F - Fernandez- Soto et al. (2004); and B - our own BPZ measurements using GOODS photometry. Entries in these two columns correspond 
to cases that do not match the (7) Adopted grism redshift, 2;grism- (8) Adopted line identification. The one case where Zspec and ^^phot are 
blank do not have good first redshift guesses, hence column (8) is left blank and the measured emission line wavelength, A in A is given 
in (7). (9) Rest frame equivalent width in A. For the one case where 2;grism is undefined, the observed EW is listed. (10) The logarithm 
of the measured line flux in erg cm~^ s"-*^. (11) Gives the methods that detected the line emission: A - aXe selection, B - blind grism 
selection. When a line is identified by both methods then the data from columns (7) - (9) are taken from method B. (12) Notes: a - two 
lines identified; b - multiple line-emitting knots; c - z^j-ism differ depending on 2;phot source; d - large 2;phot errors result in ambiguous line 
ID: e - 2;spec and ZgYism do not agree; f - no line identification allowed within 2;phot error bars, nearest expected line chosen; g - photometry 
not from GOODS release rl.lz, see individual object notes. Blank entries for columns (l)-(7) and (21) occur for the second emission line 
in sources with two emission lines (note a). These blank entries thus have the same values as for the previous line 
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TABLE 3 

Redshift, absolute magnitude, and EW statistics of ELGs split by line identification 



Property 


Ha 


[OHI] 


[Oil] 




Grism-selected LLGs (this study) 






Number 


o 




1 Q 

io 


ileasniit, z 








m mi Tm 1 m 

ililiilillLlili 


0.098 


0.308 


0.717 


25th percentile 


U. iZO 




i.UUo 


medisin 


0.132 


0.478 


1.073 


T5t]i percentile 


u.ouo 


U. O ( u 


1 '^'^7 

l.OO ( 


m n"vi rm 1 m 
iiid-A-iiii Liiii 


0.446 


0.858 


1.424 


B band absolute mag, [ABmag] 








minimum 


—20.92 


—20.97 


—23.46 


25tli percentile 


1 8 Q7 


— zu. ly 


91 R7 


median 


-18.29 


-18.93 


-21.27 


TBfVi TiPrcpnt ilp 


— 14.17 


— 18.00 


— 19.93 


maximum 


-13.03 


-14.32 


-16.61 


rest frame equivalent width, EW [A.] 








Tn mi Tni 1 m 

ililiiliilLliii 


58 


8 


22 


25th percentile 


1 no 


oy 




m pH inn 


140 


167 


107 


1 OijLl pciCclilllc 


159 


247 


198 


maximum 


223 


6418 


578 




CFRS (Hammer et al. 1997) 






Number 


95 


175 


270 


rest frame equivalent width, EW [A.] 








m mi Tm 1 Tn 

ililiilillLlili 


4 


0.7 


1.3 


ZiOiiLi peiceiiiiie 


97 


7 




median 


41 


15 


27 


75th percentile 


63 


28 


39 


maximum 


1520 


1022 


981 




SPS (Teplitz et al. 2003a) 






Number 


18 


33 


78 


rest frame equivalent width, EW [A] 








minimum 


19 


13 


6 


25th percentile 


75 


55 


45 


median 


103 


124 


68 


75th percentile 


148 


264 


117 


maximum 


394 


1479 


750 



^ There are 26 [OHI] emitters, however three are missing GOODS photometry (see Table [2]). Therefore we use the 23 sources with GOODS 
photometry to compile the Mb statistics, while all 26 sources are used to compile redshift and EW statistics. 



TABLE 4 

Color properties of HDFN field galaxies and ELGs 



field galaxies Emission-line Galaxies 



Color 
(1) 


Zphot range 
(2) 


N 
(3) 


25th 
(4) 


median 
(5) 


75th 
(6) 


Lines 
(7) 


N 
(8) 


25th 
(9) 


median 
(10) 


75th 
(11) 


^435 - i775 


< 1.62 


640 


0.63 


1.04 


1.40 


all 


45 


0.52 


0.87 


1.06 


B435 - ^775 


< 0.30 


145 


0.44 


0.70 


1.04 


Ra 


8 


0.79 


0.93 


1.11 


^435 — ^775 


0.30 - 0.76 


243 


1.07 


1.27 


1.67 


[OIII],H^ 


28 


0.91 


1.03 


1.14 


^435 — ^775 


0.76 - 1.62 


261 


0.53 


0.87 


1.31 


[OH] 


13 


0.39 


0.46 


0.73 


^775 - Z850 


< 1.62 


640 


0.07 


0.21 


0.37 


all 


45 


0.00 


0.13 


0.24 


^775 — ^850 


< 0.30 


145 


-0.03 


0.06 


0.16 


Ra 


8 


-0.02 


0.05 


0.17 


^775 — ^850 


0.30 - 0.76 


243 


0.06 


0.17 


0.27 


[OIII],H^ 


28 


0.00 


0.13 


0.20 


^775 — ^850 


0.76 - 1.62 


261 


0.24 


0.35 


0.51 


[OH] 


13 


0.13 


0.23 


0.37 



Note. — Color statistics are given for two samples of galaxies in the HDFN: field galaxies, selected purely by ^^pi^ot 

("Capakll200j), and 

ELGs selected from the grism data presented here. N (columns 3 and 8) gives the number of galaxies with GOODS photometry in each 
sample which match the Zp]^ot criteria of column (2) or contain the emission lines listed in column (7). Columns (4) and (9) give the first 
quartile color of the samples; columns (5) and (10) list the median color and columns (6) and (11) list the third quartile color. 



16 



Meurer et al. 




Grism lmag« tuts. Grfsm tatolc^ source r 287 

1200 
laaa 
Baa 



-200 

aeOO £650 iflCO £950 

Channel 

Direct Image cuts. Grfsm catalog source : 287 

2aaa 

1SGG 

I 1000 
500 
G 

aeoo £950 isco 2950 

Channel 

Fig. 1. — Steps in the processing of the grism and detection image for finding ELGs in the grism images and measuring their 
properties. Panels a and b show a 50x50 pixel cutout of the grism image before and after (respectively) subtracting a 13x3 
pixel median filtered version of the image (high-pass filtering). These cutouts are centered on an emission- line candidate selected 
from the high-pass filtered grism image. Panels c and d show cutouts of the direct image before and after high-pass filtering. 
The width of the cutout is selected to include the full range over which the direct image counterpart to the source seen in panel 
b may reside. Panel e shows the ID spectra made by extracting and summing five rows centered on the emission line from the 
grism image before (black line) and after (blue line) high-pass filtering. Panel f shows the same thing for the ID cuts through 
the direct image. The shaded regions show the pixels in the collapsed region that belong to an object found by SExtractor (the 
working object identication number is shown). 
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1 GOODSN J 12364 1.63+ 621 132.1 



2 GOODSN J 123644.75 + 62 11 57.4 



3 GOODSN J 123633. 16 + 62 1344.0 




13 ACS J1 23632.69 + 621 239.1 

0.433 mm 



16 GOODSN J1 23650.79 + 621 221. 7 



[ , 01111 , 



19 GOODSN J1 23637.76+621 235.6 




22 GOODSN J 1 23655.58+62 1 400.3 







mm 






5 GOODSN J 123648.30 + 62 1426.9 




20 ACS J1 23636.58 + 621 336.8 




0.478 




. .roiii.i 





23 GOODSN J1 23645.32 + 621 143.2 




6 GOODSN J123651.72 + 621220.5 

6.302 ■ 




24 ACS J1 23657.48 + 621 21 2.0 



Fig. 2. — Postage stamp images, 5" on a side, of the 46 ELGs, arranged by redshift. For each galaxy, two versions of the 
detection image are shown. The left panels show the detection image with spectrum extraction aperture locations containing 
emission lines marked - squares from method A, and circles from method B. The right panels shows the high-pass filtered 
detection image (see Sec 13.21 and Fig. [1]) . Here, the grism redshift is noted at the top of the panel, and the identified lines are 
noted at the bottom. The final panel shows the orientation of the images. 
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Fig. 2.— 



Continued 
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Grism sensitivity 
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Fig. 3. — Com parison of the A distribution of the hnes found in the ELGs with the grism sensitivity curve. The grism sensitivity curve 
of lWalsh Sz Pir zkal (2004) is shown in the top panel. The units of the ordinate are erg cm~^ s"-*^ A"-*^ per DN s"-*^. The middle plot shows 
the histogram of A values found with method A. The dashed lines delimit the A search range for emission lines. The bottom panel shows 
the A histogram for the method B selected ELGs. 
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(a) aXe selection 



z 10 



20 ; 
1 5 - 
10 ; 
5 - 

L_ 

0.0 ^.z , , 2.0 - 

log(Signal/NoIse) 



II 



z 10 



5 -16.0 -15.5 -15.0 

og(Flux [erg cm"^ s"^]) 

(b) blind grisnn innage selection 



1 1 
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log(Signal/Noise) 



-16.0 -15.5 -15.0 

og(Flux [erg cm~^ s~^]) 

(c) SIIS Parallel Survey 



-16.0 -15.5 -15.0 
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Fig. 4. — Line flux distribution of the ELSs found with the direct image sele ction (panel a) and grism image selection (panel b) compared 
to the ELSs found in the STIS Parallel Survey (panel c ITeplitz et al.|[2003al ). The inset in panels a,b shows the distribution of the line 
signal-to-noise ratio. 
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N = 46 

ptiles: 22.8 23.9 24.9 
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Hawaii Z3 compilation (Gowie et al. 2004) 
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28 



N = 802 

ptiles: 23.9 25.0 25.7 



I775 [ABmag] 

Fig. 5. — His togram of 2775 magnitudes of the grism-selected ELG sample (top panel) compared with the spectroscopic redshift sample 
of lCowie, et al.l (2004) (middle panel) and the photometric redshift sample of Capak (2004) (bottom panel). Objects from the latter two 
samples are only included in the histograms if they match with sources cataloged in our detection image. All 2775 magnitudes are derived 
from our ACS images. In the upper left of each corner we report the total number of sources in the sample and the 25th, 50th (median), 
and 75th percentile 2775 magnitudes. In the top panel sources with ambiguous Zgy-ism estimates (notes d-g in Table [2]) are indicated with 
the shaded histogram. 
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Fig. 6. — Rest frame B band absolute mag histograms of the the ELGs spht by hne identification, showing galaxies detected in HQ;(panel 
a); [OIII] (panel b); and [Oil] (panel c). The thick gray line shows the absolute magnitude corresponding to m = 27 ABmag for a line 
found at A = 7500A; this is a crude estimate of the faintest galaxies we are likely to find. The broken vertical lines indicate the knee of the 
luminosity function, of field galaxies at similar redshift to the ELG samples. In panel (a) the dashed line shows derived from the 
Sloan Digital Sky Survey luminosity function (Blanton, et al. 2003) while the dotted line shows the of the Two Degree Field Galaxy 
Redshift Survey (Madgwick et al. 2002). The mean 2; ~ 0.1 for both of these surveys. In panels (b) and (c) we show derived from the 
K band selected GOODS luminosity function at a mean z of 0.46 and 0.97, respectively (^Dahlen et al.|[2005r ). 
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Fig. 7. — Histogram of rest frame line equivalent widths. Histograms are normalized by the total sample size and offset vertically to ease 
comparison. In all cases our results are shown as the thick black line at top In panel (a) we compare the widths of all lines with two local 
samples: the KISS red surveys (|Salzer et al.n2001l : IGro nwall et al. 2004), shown in green (middle), and the SINGG survey (Me urer et all 
1 20061 V shown in red (bottom). The remaining panels (b,c,d) split the sou rces by lin e identification and our sample is compared to two 
surveys that extend out to mod erate redshifts: the S TIS Parallel Survey (|TeDlitz et al.ll2003al) shown with the blue middle line and the 
Canada-France Redshift Survey ([Hammer et al.l|l997l ) shown with the red line at bottom. 
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Fig. 8. — The B435 — 2775 versus 2775 — 2:850 two color plot of sources in the HDFN, using GOODS photometry (|Giavalisco et al.l 
|2004a). Small circles are sources with ^^pj^ot 

< 1.62 dCapak 2004), filled symbols are ELGs identified in this study. The symbols are color 
coded by 2:phot cind line identification respectively, where red corresponds to 2;phot < 0-30 and Ho; emitters (circles), green corresponds to 
0.3 < Zpi^ot < 0.76 and [OIII] emitters (diamonds), and blue corresponds to 0.76 < 2;phot < 1-62 and [Oil] emitters (triangles). 
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Fig. 9. — The {U — V) versus {B — z) photometry for sources in the HDFN using the Hawaii group photometry (|CaDakll2004I V The hue 
identifications are fairly well sorted in this plane. The symbols used are the same as in Fig. [8] 
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Fig. 10. — Comparison of spectroscopic redshifts from fCowie. et al.|[2'004l ) and grism redshifts from this work. Measurements from 
method A are shown with sohd symbols, measurements from method B are shown as open symbols. The symbol shape and color indicate 
the grism line identification: Ho; emitters are red circles, [OIII] emitters are green triangles, and [Oil] emitters are blue squares. The unity 
relationship is shown as a solid line, sources outside the dashed lines at = ±0.105 are considered outliers. 
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Fig. 11. — Angular size plotted against redshift for the sources in our field. Sources with photometric redshifts are plotted as small black 
dots. Yellow filled circles mark sources with spectroscopic redshifts. The large color filled symbols mark the sources with grism redshifts; 
red filled circles indicate Ha emitters, green diamonds indicate [OIII] or H/3 emitters, and blue triangles indicate [Oil] emitters. 



